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The quantum anomalous Hall effect (QAHE) is a fundamental quantum transport phenomenon 
that manifests as a quantized transverse conductance in response to a longitudinally applied electric 
field in the absence of an external magnetic field, and promises to have immense application po¬ 
tentials in future dissipation-less quantum electronics. Here we present a novel kinetic pathway to 
realize the QAHE at high temperatures by n-p codoping of three-dimensional topological insulators. 

We provide proof-of-principle numerical demonstration of this approach using vanadium-iodine (V-I) 
codoped Sb 2 Te 3 and demonstrate that, strikingly, even at low concentrations of ~2% V and ~1% I, 
the system exhibits a quantized Hall conductance, the tell-tale hallmark of QAHE, at temperatures 
of at least ~ 50 Kelvin, which is three orders of magnitude higher than the typical temperatures at 
which it has been realized so far. The proposed approach is conceptually general and may shed new 
light in experimental realization of high-temperature QAHE. 

PACS numbers: 73.20.At; 74.43.Cd; 75.70.Tj. 


Two-dimensional electron systems (2DESs) are versa¬ 
tile playgrounds for complex quantum transport phenom¬ 
ena that are of interest both from fundamental and ap¬ 
plication points of view. One prominent example is the 
quantum Hall effect (QHE) [lj, [2|, a quantum analogue 
of the classical Hall effect, whereby the application of 
a magnetic field to a 2DES results in the quantization 
of the transverse conductance and a vanishing longitudi¬ 
nal conductance. Crucially, the mechanisms underlying 
this phenomenon are the formation of an insulating bulk 
and remarkable chiral conducting edge states (Fig. la). 
This striking unidirectionality makes backscattering im¬ 
possible, and renders transport along these edge states 
remarkably robust against any impurities |,@]. This 
robustness essentially renders the edges of 2DESs mani¬ 
festing the QHE as perfectly conducting one-dimensional 
(ID) wires, which would be of great interest as potential 
building blocks, such as interconnects between chips, in 
dissipation-less electronics. However, a major constraint 
for practical use of this fundamental quantum transport 
phenomenon is the requirement of very strong magnetic 
fields, which is impractical for realistic applications. 

Realizing the quantum anomalous Hall effect (QAHE) 
would circumvent this problem because it manifests the 
same hallmark features as the QHE, but does not require 
the application of an external magnetic field jH]. The rea¬ 
son why the QAHE does not require the application of 


an external field is because the latter’s role is played by 
an intricate cooperation between the intrinsic magnetism 
that breaks the time-reversal symmetry and spin-orbit 
coupling in the 2D (quasi-2D) insulating host material. 
Initial proposals for realizing the QAHE were based on 
honeycomb lattice models 0. Since then, especially af¬ 
ter the experimental discoveries of graphene (with an in¬ 
herent honeycomb lattice structure) [6| and topological 
insulators (TIs) [3, [§J, much effort has been made to ex¬ 
ploring new platforms for realizing QAHE in these and 
related systems IM- TIs are inherently superior to 
graphene for this purpose because of the formers’ intrinsi¬ 
cally strong spin-orbit coupling, which narrows the search 
for suitable materials to those that could harbor intrin¬ 
sic ferromagnetism. One possible way to introduce fer¬ 
romagnetism in TIs is by introducing magnetic dopants 
1121 l2l(|23j. similar to the approach taken in realizing di¬ 
luted magnetic semiconductors [24]. This approach has 
led to the first experimentally realized QAHE systems in 
Cr-doped (Bi,Sb) 2 Te 3 thin films [25l427t ] . Yet to date, all 
the experimentally observed QAHE were measured only 
at extremely low temperatures, which would still pro¬ 
hibit potential applications of this fundamental quantum 
transport phenomenon. 

In this Letter, we propose a versatile approach to re¬ 
alize high-temperature QAHE by using compensated n- 
p codoping to establish high-temperature ferromagnetic 
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FIG. 1: (a) Schematic of a Hall-bar device with chirally prop¬ 
agating edge modes, (b) The 4x4x1 supercell of Sb 2 Te 3 with 
only one of three QLs shown for clarity. The notations are: 
VI is the site, where one V atom substitutes one Sb atom; 
V2-A - V2-D are sites that can be occupied by a second V 
dopant; I-A - I-D designate sites that can be occupied by one 
I dopant, (c) V-V interaction energy E re \ vs. the V-V sepa¬ 
ration dy-v of V doped (open black circles) and V-I codoped 
Sb 2 Te 3 (filled red squares), respectively. For the V doped 
case, the label of each open circle indicates the position of the 
second V atom. In the V-I codoped case, one I is fixed at I-A, 
and the two labels in each square bracket denote the positions 
of the second V and second I, respectively, (d) Magnetic cou¬ 
pling between two V dopants vs. dv-v of V doped (open 
black circles) and V-I codoped (filled red squares) Sb 2 Te 3 , 
respectively. 


order in 3D-TIs. This doping approach in large part 
does not alter the bulk gaps and thus addresses the pri¬ 
mary challenge that has hindered experimental realiza¬ 
tion of high-temperature QAHE in TIs so far 2^-27|. 
We provide numerical proof-of-principle demonstration 
of this approach using a prototypical example, by codop¬ 
ing Sb 2 Te 3 with magnetic p-type vanadium (V) and non¬ 
magnetic n-type iodine (I) dopants. These dopants are 
chosen because of their preference to form n-p pairs due 
to their mutual electrostatic attraction, thereby enhanc¬ 
ing the solubility of both dopants. More importantly, 
we show that whereas doping with V alone would unde¬ 
sirably shift the Fermi level into the valence band and 
shrink the bulk band gap, codoping with I recovers the 
insulating nature while also significantly restores the in¬ 
trinsic bulk band gap of the TI. Moreover, for all the 
codoping concentrations considered in our study, the cou¬ 
pling between the hard magnetic V dopants always gives 
rise to a stable ferromagnetic order. Strikingly, even at 
a low codoping concentration of 2.08% V and 1.39% I 
(which amounts to compensated one-to-one V-I codop¬ 
ing), an estimated Curie temperature about 49 Kelvin 


FIG. 2: Band structures of V doped and V-I codoped Sb 2 Te 3 . 
(a) Pure Sb 2 Te 3 with a bulk band gap of 123 meV; (b) For 
2.08% V doping, the Fermi level shifts into the valence band, 
and the bulk gap decreases to 76 meV; (c) Additional 0.69% 
I codoping slightly moves the Fermi level towards the bulk 
gap and enlarges the gap to 81 meV; (d) A further increase 
of I codoping to 1.39% compensates the V dopants and shifts 
the Fermi level into the bulk gap of 89 meV, effectively recov¬ 
ering the insulating nature of the system. (d-l)-(d-4) Differ¬ 
ent characters of the bands shown in panel (d), obtained by 
projecting the Kohn-Sham states onto the local orbitals of a 
single atom for each element. 


can be reached. Aside from the bulk properties, the band 
gaps associated with the surface states of the correspond¬ 
ing thin film structures can also be as large as 53 meV, 
much larger than what is needed to sustain ferromagnetic 
order at the high Curie temperatures. Explicit calcula¬ 
tions of the Berry curvatures further confirm that the 
V-I codoped Sb 2 Te 3 thin films can harbor QAHE with 
the Hall conductance quantized to be a xy = e 2 /h . The 
proposed approach thus offers a versatile route to the re¬ 
alization of high-temperature QAHE in different TIs and 
related materials. 

Our first-principles calculations were performed us¬ 
ing the projected augmented-wave method [28( as im¬ 
plemented in the Vienna Ab-initio Simulation Package 
(VASP) [29[. The generalized gradient approximation 
(GCA) of Perdew-Burke-Ernzerhof (PBE) type was used 
to treat the exchange-correlation interaction [3C|. A 
4x4x1 Sb 2 Te 3 supercell is chosen to study the magnetic 
coupling between a pair of V dopants. Each supercell is 
composed of three quintuple layers (QLs), including 96 
Sb atoms and 144 Te atoms. For thin film calculations, 
the film thickness was chosen to be six QLs, including 
192 Sb and 288 Te atoms. A vacuum buffer space of 
30 Awas used to prevent the coupling between adjacent 
slabs. The kinetic energy cutoff was set to 250 eV. Dur- 
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ing structural relaxation, all atoms were allowed to re¬ 
lax until the Hellmann-Feynman force on each atom is 
smaller than 0.01 eV/A. The Brillouin-zone integration 
was carried out by respectively using 3x3x2 and lxl 
Monkhorst-Pack grids for bulk and thin film systems. 
Unless mentioned otherwise, spin-orbit coupling was con¬ 
sidered in all calculations, and the GGA+U method was 
used with on-site repulsion parameter U=3.00 eV [3lj | 
and exchange parameter J =0.87 eV. The Curie temper¬ 
ature Tq was estimated within the mean-field approxima¬ 
tion ksTc = 32], where ks is the Boltz¬ 

mann constant, x is the dopant concentration, and Joi 
is the on-site exchange parameter obtained from the to¬ 
tal energy difference between the ferromagnetic and anti¬ 
ferromagnetic configurations. 

To elucidate the role of codoping in realizing high- 
temperature QAHE, we first investigate the energetics 
and magnetic behavior of Sb 2 Te 3 solely dop ed with V, 
whereby V dopants substitute Sb atoms [33[. The inter¬ 
action energy between two V atoms is defined as 


Ere 1 = E2V + Eqv — 2Eiy, 


( 1 ) 


where Ejy is the total energy of the Sb 2 Te 3 systems con¬ 
taining i (i= 0, 1, 2) V atoms. We found that the V- 
V interaction is repulsive (Fig. lc), with E re l > 0 for 
all the positions considered. Thus the energetics dictate 
that the V dopants arrange according to a diluted distri¬ 
bution. We also found that the V-V magnetic coupling 
is ferromagnetic for all the V-V separation considered in 
our study, strongly indicating a preponderance toward 
a diluted ferromagnetic (Fig. Id) - the magnetic cou¬ 
pling between two V moments is defined as the total 
energy difference between the ferromagnetic and anti¬ 
ferromagnetic configurations. 

One important finding from our comparison of the 
electronic band structures of pristine Sb 2 Te 3 (Fig. 2a) 
and V-doped Sb 2 Te 3 (Fig. 2b) is that solely doping the 
TI with V-dopants drives the system from an insulator 
(Fig. 2a) to ap-type semiconductor (Fig. 2b), even at a 
low concentration of 2.08%. The Fermi level is thereby 
shifted to the valence band (Fig. 2b), which is consis¬ 
tent with previous experimental findings [33} ■ V doping 
also decreases the bulk gap from 123 meV to 76 meV at 
the 2.08% V concentration, which originates from the V 
contribution to the valence band maximum (VBM) (Fig. 
2(d-l)) jUj]. Consequently, V doping alone cannot con¬ 
vert Sb 2 Te 3 into a magnetic TI. 

We now demonstrate that to recover the insulating na¬ 
ture of Sb 2 Te 3 it is necessary to compensate the extra 
charge introduced by V dopants, by adding an additional 
type of dopant with opposite charge, i.e., in this case an 
n-type dopant. For the purpose of our proof-of-principle 
demonstration we chose I as the dopant to substitute 
Te in Sb 2 Te 3 , which results in an n-type doping. This 
choice was inspired by the fact that SbTel and BiTel 


are naturally existing materials 135} . Moreover, I pos¬ 
sesses larger spin-orbit coupling than Te, which is likely 
to have the added benefit for enhancing the topological 
phases fdtil. [ 37 }. Using an n-p codoping approach has also 
the added benefits of stabilizing the constituent dopants 
(because of their mutual electrostatic attraction, which 
leads to the codopants forming closely spaced pairs dur¬ 
ing the doping process), increasing their solubility, and 
potentially leading to a higher magnetic transition tem¬ 
perature, as in the case of diluted magnetic semiconduc¬ 
tors 38]. Additionally, we note that the experimentally 
fabricated Sb 2 Te 3 is usually a light-doped TI (therefore 
a conductor) due to the defects, which, based on our 
compensation principle, can be potentially tuned to be 
an insulator by properly adjusting the ratio of n- and 
p-doping. 

Codoping Sb 2 Te 3 with V and I gives rise to striking 
changes in the electronic band structure. First, codop¬ 
ing with 0.69% I, in addition to 2.08% V doping, shifts 
the Fermi level upwards (Fig. 2c). For this low I-doping 
concentration, the system is still a p-type semiconduc¬ 
tor. When the I concentration further increases to 1.39% 
(which entails one-to-one V-I codoping), the system re¬ 
covers its insulating phase, and the bulk gap also in¬ 
creases to 89 meV (Fig. 2d). Further analysis of the 
characters of different elements in the electronic band 
structure (Figs. 2d(l)- 2d(4)) reveals that the valence 
band close to the Fermi level is mainly dominated by 
V, Sb, and Te atoms, while I contributes mainly to the 
lower part of the valence bands. These results directly 
reflect the strong bonding interaction between I and the 
surrounding V/Sb atoms, because I has larger electroneg¬ 
ativity than Te, and the resulting large energy difference 
between the bonding and anti-bonding orbitals helps to 
enlarge the bulk gap. Therefore, the introduction of I 
not only shifts the Fermi level into the bulk gap but also 
enlarges the gap size, which are the salient features of 
the compensated n-p codoping scheme. Taken together, 
these results demonstrate that it is possible to preserve 
an insulating state for bulk Sb 2 Te 3 even in the pres¬ 
ence of doping, provided that one-to-one compensated 
V-I codoping is used. 

Preserving the insulating phase during doping is one 
key requirement for realizing high-temperature QAHE. 
The second requirement is to achieve a strong ferromag¬ 
netism. We will therefore examine the magnetic prop¬ 
erties of the V-I codoped system, and in particular elu¬ 
cidate under what conditions ferromagnetic order is fa¬ 
vored. Before doing so, it is important to investigate in 
more detail how two separate V-I pairs interact in the 
host system. Our calculations show, first, that the pair- 
pair interaction energies E ve 1 are attractive (Fig. lc). A 
second conclusion from our calculations is that the mag¬ 
netic coupling between two magnetic moments associated 
with V dopants is ferromagnetic, with AUfm-afm < 0, 
for all the considered dopant separation ranges. In ad- 
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dition, we found that, when compared to the pure V- 
doped case, the resultant magnetic moment of a V-I pair 
is larger, but the coupling strength between two such 
pairs is somewhat smaller. This can be attributed to 
the fact that the V-I codoped system becomes insulat¬ 
ing, whereas the V-doped system is metallic and has de¬ 
localized charge carriers that mediate stronger magnetic 
coupling. Taken together, these findings indicate that 
the n-p codoping scheme efficiently leads to a stable fer¬ 
romagnetic state even at relatively low codoping concen¬ 
trations, which is highly advantageous compared to other 
approaches that rely on heavy magnetic doping [25l - l27j . 

The ferromagnetism induced using our codoping 
scheme is highly robust and has the potential to exist 
at high temperatures. For the particular doping concen¬ 
trations in our proof-of-principle demonstrations, using 
the magnetic coupling strengths extracted from our first- 
principles calculations, we estimate the Curie tempera¬ 
ture within mean-field theory [32] to be Tq = 49 (152) 
Kelvin for codoping concentrations of 2.08% V and 1.39% 
I, respectively [the results were obtained within GGA+U 
(GGA) by placing two V-I pairs in diluted configura¬ 
tion at one QL], With further increase of the codoping 
concentration it is possible to raise the Curie tempera¬ 
ture, as long as the dopant pair concentrations are still 
low enough to be away from potential clustering of the 
dopants [391 ]. Therefore, the present estimations give a 
lower limit on the Curie temperature at such a low codop¬ 
ing concentration. 

Our proof-of-principle demonstration thus shows that 
V-I codoped Sb 2 Te 3 is a ferromagnetic insulator with a 
sizable gap, making the system a highly desirable can¬ 
didate for realizing high-temperature QAHE. To demon¬ 
strate that such systems can indeed harbor QAHE for 
a wide range of codoping concentrations, we have calcu¬ 
lated and analyzed the band structures of intrinsic and 
V-I codoped Sb 2 Te 3 thin films and shown that these sys¬ 
tems yield the quantized Hall conductance a xy = e 2 //i 
expected for the QAHE, i.e., without the external mag¬ 
netic held. To that end, first, we have shown that these 
codoped systems harbor large energy gaps around the 
Dirac point. In the absence of dopants, the surface of a 
Sb 2 Te 3 thin him hosts a massless Dirac fermion, mani¬ 
fested as the trademark linear Dirac dispersion near the 
gamma point [itj (Fig. 3a). After introducing magnetic 
dopants, the Dirac fermion acquires a hnite mass due 
to the presence of the intrinsic ferromagnetism, which in 
principle can give rise to the QAHE j§]. Here we show the 
calculated band structures and Berry curvatures, along 
the high-symmetry directions near the gannna point, for 
three different V-I codoping concentrations 2.08% V and 
1.39% I (Fig. 3b), 4.16% V and 2.78% I (Fig. 3c), and 
6.25% V and 4.17% I (Fig. 3d). For all the codoping 
concentrations the Fermi level lies inside the gap, which 
is enlarged from 53 meV (Fig. 3b) to 84 meV (Fig. 3d) 
when the codoping concentration increases. This proves 
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FIG. 3: (a) Band structure of a Sb 2 Te 3 thin Him calculated in 
a 4 x 4 Sb 2 Te 3 supercell with a thickness of 6QLs along high 
symmetry directions, (b)-(d) Solid lines: Band structures 
of V-I codoped Sb 2 Te 3 with the film thickness of 6QLs at 
different codoping concentrations: (b) 2.08% V and 1.39% 
I, (c) 4.16% V and 2.78% I, and (d) 6.25% V and 4.17% I. 
Square symbols: Corresponding Berry curvatures. 


that the V-I codoping drives the conducting surface into 
an insulator with an appreciable surface band gap, allow¬ 
ing the formation of high-temperature QAHE. 

A separate and unambiguous evidence for the mani¬ 
festation of QAHE is obtained by integrating the Berry 
curvature of the occupied valence bands using the expres¬ 
sion 


4l|, [421 


n(k) 


53 f n 53 

n n'^n 


2Im(V’nfc | V x | t/V k ) ( i>n> k | V y | Ipnk ) 
{E n , - E n )2 


_ ( 2 ) 

where n is the band index, E n and ijj n k are the eigen¬ 
value and eigenstate of the band n, v x = dE/dk x and 
v y = dE/dky are velocity operators along x and y di¬ 
rections within the film plane, and f n = 1 for all oc¬ 
cupied bands. This integration gives the corresponding 
Hall conductance (42|. The Berry curvature distribution 
along the high symmetry directions shows a large nega¬ 
tive peak near the gamma point and zero elsewhere (Figs. 
3(b)-3(d)). As a consequence, the total integration or 
the Hall conductance with the Fermi level lying inside 
the band gap must be nonzero. With the increase of the 
surface band gap, the region with finite Berry curvatures 
broadens, and the curvature peak gradually decreases, 
potentially signifying the unchanged Hall conductance. 
Our quantitative calculations reveal that the Hall con¬ 
ductance is indeed quantized, i.e. a xy = e 2 /h, for all the 
three codoping concentrations, and thus that the system 
is in a QAHE state. 

In conclusion, our findings demonstrate that codoping 
TIs provides a versatile route to realize the QAHE at 
high temperatures (43|. Our proof-of-principle demon¬ 
strations revealed that codoping results in: 1) the for¬ 
mation of large band gaps, 2) high Curie temperatures, 
3) robust quantized Hall conductance. These manifes¬ 
tations originate from the realization of hard ferromag¬ 
netism (in our case due to V-dopants), strong spin-orbit 
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coupling (here because of I-dopants), and, most crucially, 
the pinning of the Fermi level inside the largely preserved 
(upon codoping) intrinsic TI bulk gap. Even at very low 
concentrations, we have shown that the lower-bound tem¬ 
perature for observing the QAHE in our codoped TI is 
~ 50 Kelvin, which is three orders of magnitude higher 
than what was previously observed |25H27], and by fur¬ 
ther increasing and tailoring the codoping and TI envi¬ 
ronment it is feasible to further increase this temperature. 
So far, the QAHE was experimentally observed by us¬ 
ing the host TI (Bi,Sb) 2 Te 3 doped with chromium atoms 


25-27]. These experiments could be readily extended to 


using the proposed codoping approach to realize high- 
temperature QAHE. n-p codoping, while never used in 
the context of the QAHE, has been previously used, e.g., 
by ourselves in the studies of different kinds of functional 
materials (dsl. fbl} . 
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